Non-negative Wigner functions for orbital angular momentum states 
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The Wigner function of a pure continuous-variable quantum state is non-negative if and only if the state is 
Gaussian. Here we show that for the canonical pair angle and angular momentum, the only pure states with non- 
negative Wigner functions are the eigenstates of the angular momentum. Some implications of this surprising 
result are discussed. 
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For continuous variables, the Wigner function yj] is a very 
useful tool that establishes a one-to-one correspondence be- 
tween quantum states and joint quasiprobability distributions 
of canonically conjugate variables in phase space (position 
and momentum, in the standard case). However, it can take 
on negative values, a property that distinguishes it from a true 
probability distribution yj, |3l [4]. Indeed, this negative char- 
acter is associated with the existence of quantum interference, 
which itself may be identified as a signal of nonclassical be- 
havior H. 

In consequence, the characterization of quantum states that 
are classical, in the sense of giving rise to non-negative 
Wigner functions, is a topic of undoubted interest. Among 
pure states, it was proven in a classicalpaper by Hudson ]6[] 
(later generalized by Soto and Claverie H to multipartite sys- 
tems) that the only states that have non-negative Wigner func- 
tions are Gaussian states 18t|9[]. This is one of the main reasons 
for the prominent role these states play in modern quantum in- 
formation tfToh - 

The original definition of the Wigner function has also been 
extended to discrete systems (see Ref. ifTTIl for a comprehen- 
sive review). Again, the classification of states with non- 
negative Wigner functions is an amazing problem that has 
been solved quite recently by Paz and coworkers lfl2i [l3ll and 
Gross 03. H, so that the role of Gaussian states is now taken 
on by stabilizer states. Interestingly, these are the only states 
that can be simulated efficiently in classical computers iflcill . 

Between these two cases (whose proofs are otherwise com- 
pletely different), we have the interesting situation of canon- 
ical pairs, such as the angle and orbital angular momentum 
(OAM), for which one variable is continuous while the other 
one is discrete lfl7ll . The associated phase space is the discrete 
cylinder S\ x Z, where S\ stands for the unit circle (associ- 
ated to the angle) and the integers Z translate the discreteness 
of the OAM. The physical example we have in mind is the 
OAM of photons. This is an emerging field that has given 
rise to many developments, ranging from optical tweezers to 
high-dimensional quantum entanglement, or fundamental pro- 
cesses in Bose-Einstein condensates, to cite only a few rele- 
vant examples lfl8ll . 

The seminal paper of Allen et al. lfl9ll firmly established 
that the Laguerre-Gauss modes carry a well-defined OAM. 
They appear as annular rings with a zero on-axis intensity 
and an azimuthal dependence cxp(i^0) that gives rise to spi- 



ral wave fronts. The index £ takes only integer values and can 
be seen as the eigenvalue of the OAM operator. Since then, 
several methods have been established to produce light beams 
with the required azimuthal phase structure, among these spi- 
ral phase plates, forked holograms, and spatial light modula- 
tors are perhaps the most versatile. In this way, a variety of 
modes with helical phase fronts but different transverse pat- 
terns (such as Bessel, Mathieu, or hypergeometric beams) can 
be routinely generated in the laboratory [20;]. 

The goal of this work is precisely to determine the pure 
states of these OAM-carrying systems for which the Wigner 
function is non-negative, filling in this way a long overdue 

gap- 
To be as self-contained as possible, we first introduce some 
basic notions for the problem at hand of cylindrical symmetry. 
We are concerned with the planar rotations by an angle <j> gen- 
erated by the angular momentum along the z axis, which for 
simplicity will be denoted henceforth as L. We do not want 
to enter in a long discussion about the possible existence of 
an angle operator ll2lll . For our purposes here, the simplest 
solution is to adopt two periodic angular coordinates, e.g., co- 
sine and sine, that we shall denote by C and S to make no 
further assumptions about the angle itself. One can concisely 
condense all this information using the complex exponential 
of the angle E = C + iS, which satisfies the commutation 
relation 



[E, L]=E. 



(1) 



In mathematical terms, this defines the Lie algebra of the 
two-dimensional Euclidean group E(2), which is precisely the 
canonical symmetry group for the cylinder. 

The action of E on the basis of eigenstates of L is E\l) = 
\£ — 1), and it possesses then a simple implementation by 
means o f a p hase mask removing a charge +1 from a vor- 
tex state I22ll23ll . Since the integer I runs from —00 to +00, 
E is a unitary operator whose eigenvectors 



1 



(2) 



form a complete basis and describe states with well-defined 
angle. In the representation generated by them, L acts as — id,/, 
(in units of K = 1). 
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Given the key role played by the displacement operators in 
settling the Wigner function for the harmonic oscillator, we 
introduce a unitary displacement operator 



D(£,</>) 



E-e e -*4>L 



(3) 



where a(£, <fi) is a phase required to avoid plugging in extra 
factors when acting with D. The conditions of unitarity and 
periodicity restrict the possible values of a, although a sen- 
sible choice is a(£, 4>) = —£<fi/2. Note that here we cannot 
rewrite Eq. (01 as an entangled exponential, since the action 
of the operator to be exponentiated would not be well defined. 

We use as a guide the analogy with the continuous case and 
introduce the mapping rf24ll 



W s (l,, 



Tv[gw( 



(4) 



which maps the density operator into a Wigner function via 
a kernel w defined as a double Fourier transform of the dis- 
placement operator l25ll : 

^ ' fez 1 ' 2 *' 

(5) 

where the integral extends to the 2ir interval within which the 
angle is defined. This mapping is invertible, so one can recon- 
struct the density operator as 



g = 2tt 



E 



Ws 



(6) 



2jt 



The (Hermitian) Wigner kernels w(£, <j>) are a complete or- 
thonormal basis (in the trace sense) for the operators acting on 
the Hilbert space of the system. In addition, they are explicitly 
covariant; i.e., they transform properly under displacements, 
w(£, 4>) = £>(£, <ft) w(0, 0) D*(£, 4>). In fact, these properties 
guarantee that the Wigner function defined in Eq.© bears all 
the good properties required for a probabilistic description. In 
particular, it reproduces the proper marginal distributions, that 
is, 



ev ., V) - \ Yl Q\(p), I W s (£, 4>)d<j) = (£\g\£) 

2tt 



(7) 

Finally, the overlap of two density operators is proportional to 
the integral of the associated Wigner functions: 



Tr^cx^T I W & 



(e,<t>)w & (i,(i>) dct>. 



(8) 



This property (often called traciality) offers practical advan- 
tages, since it allows one to predict the statistics of any out- 
come, once the Wigner function of the measured state is 
known. 

We remark that this approach to the Wigner func- 
tion is grounded in the axiomatic method developed by 
Stratonovich JH] and Berezin JH (see also Ref. J27tl). It 
is possible to follow alternative routes, such as, introducing 



a Wigner function as the Fourier transform of some general- 
ized characteristic function l28ll . This has been pursued also 
for the group E(2) rf29tl . However, these apparently disjoint 
formulations turn out to be equivalent for most practical pur- 
poses i30ll . 

To give an explicit form of the Wigner function © we need 
to evaluate it in a basis. Using the OAM eigenstates, we get 



re: 



2tt 2 4-^ £" + 1 



/2 C 



-(2£' + l)i<p 



t',l"& 

x (£ + £" -£'\g\£ + £" + £' + 1) . (9) 

This looks rather cumbersome due to the second sum in 
Eq. © and sometimes is preferable to work in the angle rep- 
resentation, for which one easily finds 



W- e {£,ct>) 



1 

2^ 



-4>' /2\ Q \4> + cp' /2)e^' i d<t>' . (10) 



This coincides with the result of Mukunda 13 ll 13211 (see also 
Ref. rf33ll ) and bears a resemblance with the standard Wigner 
function for position and momentum that is more than evi- 
dent. Note that using this latter function in terms of transverse 
coordinates, as is often done in classical optics l34ll . is not ap- 
propriate for the geometry of the cylinder, which is the natural 
domain in which the Wigner function should be defined. 

We have now all the ingredients needed to accomplish 
our program. In what follows, the Fourier transform of 2tt- 
periodic functions (i.e., with domain in Si), defined as 



^ 9 ){k) = ^- j <?(</>) 



(11) 



with k G Z, will play a relevant role. We first state our main 
result, which can be viewed as analogous to the Hudson theo- 
rem for the canonical pair angle and angular momentum. 

Theorem (Classical OAM states). The Wigner function of a 
pure state \tp) is non-negative if and only if \tp) is an OAM 
eigenstate \£q}. 

Proof. The sufficiency is obvious since the Wigner function 
for the state \£ Q ) is W\e )(£, 4>) = 5u /(2n). The delicate 
point is to prove the necessity. Before proceeding, we sketch 
the idea behind the proof. The first step is to show that the 
wave function [and thus, the integrand in Eq. (fTQb l must be of 
constant modulus. The second step is then to corroborate that 
the Wigner function can only be non-zero for a single value 
of £. Traciality permits us to derive an equation that shows 
that this value of £ cannot vary over <j>, and that indeed the 
only states with non-negative Wigner functions are the OAM 
eigenstates. We start with the following lemma. 

Lemma 1. If the Fourier transform of a smooth, complex, 2tt- 
periodic function g(4>) is non-negative, then the integration 
kernel g(cj> — (f>') is non-negative. 
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Proof. By a direct calculation we can check that 

/ ff(^-0 / )e-^' fc # , =2 7 r(^)(fe)e- i ^ ) (12) 

J2-IX 

so, for any smooth test function x(4>) = J2kei X(k) e ~ > it 
holds 

/ X *W>) g{4>-4>') x{4>') <W = 4^ 2 J2 l^( fc )l 2 (^)(*) ■ 

1/2,7 fcez 

(13) 

It is clear that the non-negativity of the kernel g(<j> — cf>') 
follows from the non-negativity of the Fourier transform 
(Tg){k). □ 

We apply the lemma to 

x(<t>) = \[&**(<l> - ci) + M0 - Pa)] , (14) 

Here, (5 27 r denotes the periodic delta function (or Dirac comb) 
of period 2-k and c\ , c 2 G 5i . For this function we have 
lx(^)| 2 = {1 + cos[fc(ci - c 2 )]}/(87r 2 ), so the sum in the 
right-hand side of Eq. $13[ reduces to 

g(0)/2 + [g( Cl - c 2 ) + g(c 2 - d)]/4 . (15) 

Consequently, for a function g(<j!>) whose Fourier transform is 
non-negative, the kernel g(<f> — <fr') must also be non-negative 
on the test functions ( TPfl i for all the possible parameters 

ci,c 2 6 Si. 

For a pure state \ip), the Wigner function (fTOb is just the 
Fourier transform of ip*(</> + <j)' /2) ip(4> — </>'/2), where we 
have expressed the wave functions in the angle representation. 
By Lemma 1, for the test functions (fl4l i the non-negativity of 
W\m leads to 

|^)| 2 > 1^(0-0/2)11^(0 + 0/2)1, (16) 

with a = c\ — c 2 . This implies that 1^(0)1 cannot have any 
minima and the modulus of ip must thus be flat over Si. 
To proceed further we need a technical detail. 

Lemma 2. If a function f(k) : Z — > C has an inverse Fourier 
transform of constant modulus over <p, then 

fcez 

Proof. Let us first introduce the operator 

A = E f(m-k)\m)(k\. (18) 

m.fcGZ 

One can check that it can be expressed in a diagonal form in 
the angle basis, namely 

A= [ O("0) (19) 

If (J r_1 /)(0) has constant modulus, it can be written as 
= c e iX ^\ where A is a real function. Therefore, 



we have A A* = \c\ 2 1. But according to the definition ( fT8l . 
this is tantamount to the orthogonality relation 

(■n\m)(k\f(m-k)\k')(m'\f*(m'-k')\n+j) =0. 

m,fc£Z m' ,k' £Z 

The Plancherel formula allows one to cancel the diagonal 
parts, so we are led to 

£/(n-fc)/*(« + J-fc) = 0, (20) 

fcez 

whence the result follows. □ 

Next, for every <fi, we consider the Wigner function of the 
state as a function exclusively of the discrete index I; that is, 
f^{£) = W w (£, 0) : Z -> R (in fact, W is real valued), and 
make use of the fact that the (inverse) Fourier transform of 
ftj>{l) has a constant modulus over <f>. Then, by Lemma 2, the 
orthogonality 

E/*W/;(^ + O=0. W'/O, (21) 

^ez 

must hold for all <f> G S\. But since / is non-negative on the 
whole phase-space, this is only possible if / is equal to zero 
for all but one £q. Note that, in principle, £q may depend on <f>. 
Taking into account the marginal distribution Q, we see that 
W(£,(b) = 5 UoW /(2n). 

We now make use of the fact that the state \tp) is pure [that 
is, Tr(p 2 ) = 1]. From the traciality property, one can show 
that the Wigner function representing the product of two den- 
sity operators g and a can be expressed as 

W 6 »{£,4>) = ±- I W^ + £i,0 + Vi/2) 

^ « 1 £ 2 GZ"' 27r 

x Wa(£ + 12,0 + V2/2) e l ^~ 1 ^ #i# 2 • (22) 

We apply this to the pure state \ip) whose Wigner function is 
oftheform<5 m)(0) /(27r). 

Without loss of generality, we can assume that Iq ((f> = 0) = 
and may revert this choice later by a displacement — > 
D(£ , 0)|0)- Then, Eq. <|22) becomes 

= — L-/ e^^/^-^/^d^d^. (23) 

This means that the integral of the imaginary part must vanish, 
while the integral of the real part must be equal (27r) 2 . This 
is only possible if the exponential is exactly one for all the 
arguments (-01, ^2); i.e., £a(ipi/2)ip2 = £0(^2/2) ipi mod 
2ir. This is only possible when Iq = 0. □ 

We have shown that if the Wigner function of a pure state is 
non-negative, then it is necessarily a Kronecker delta, and thus 
stems from an OAM eigenstate, which concludes the long yet 
instructive proof of our theorem. 
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It is worth stressing that for the continuous case the no- 
tions of coherent states, Gaussian wave packets, and states 
with non-negative Wigner functions (often identified as non- 
classical states) are completely equivalent. However, special 
care must be paid in extending these ideas to other physical 
systems like OAM, since they lose their equivalence. 

For example, OAM coherent states \£o,<fio) m the cylin- 
der ll35ll can be expressed in the angle representation by 



Wo, 



D ifo(0-0o) 



^3 (Oil) 



1 



where $3 denotes the third Jacobi theta function. However, 
despite the key role played by this function in angular prob- 
lems, a simple calculation Oal immediately reveals that the 
Wigner function for them takes negative values. 
In the same vein, the states 



1 



exp(Atcos< 



(24) 



whose associated probability distribution is precisely the von 
Mises distribution |2lll . are usually taken as Gaussians for this 
problem. One can easily check that their Wigner function also 
takes negative values. 



Even with all these cautions, the characterization we have 
presented of OAM eigenstates as the only ones with non- 
negative Wigner function has interest in its own, although, 
unfortunately, they cannot be viewed as Gaussian states. 

A topic of interest is the characterization of unitaries that 
preserve the non-negativity. Obviously, all the displacement 
operators are of this kind. But the exponential of an arbitrary 
real function f(L) also preserves non-negativity and this in- 
cludes quadratic exponentials, which are essential for a full 
quantum reconstruction of vortex states l36Tl . 

Finally, let us mention that a question that naturally arises is 
whether our result can be extended to mixed states. Although 
this question has been approached by using the notion of the 
Wigner spectrum [3711 and explored quite recently for contin- 
uous variables 13811 . in our case a simple extension seems dif- 
ficult and will be the object of our future work. 
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